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ABSTRACT 

We present new results from a systematic study of absorption line strengths of galaxies in clusters approaching 
redshifts of unity. In this paper, we specifically compare the strengths of the high-order Balmer absorption features 
of H7 and HS in E/SOs in the four clusters Abell 2256 (z=0.06), CLI 358+62 (z = 0.33), MS2053-04 (z = 0.58), and 
MS 1054-03 (z = 0.83). By comparing the correlation of Balmer line strength with velocity dispersion for E/SOs 
in the three clusters, we find moderate evolution in the zero-point of the (H(5/i+H7A)-cr relation with redshift. The 
trend is consistent with passive evolution of old stellar populations. Under the assumption that the samples can be 
compared directly, we use single-burst stellar population synthesis models to constrain the last major occurrences 
of star-formation in the observed E/SOs to he Zf > 2.5 (95% confidence). We have compared the evolution of 
the Balmer absorption with the evolution of the B-band fundamental plane and find that simple stellar population 
models agree very well with the data. While the best agreement occurs with a low value for ri„,, the data provide 
strong confirmation of the time-evolution in recent stellar population models. 

Subject headings: galaxies: evolution, galaxies: elliptical and lenticular, galaxies: clusters: individual (Abell 
2256,CL1358+62,MS2053-04,MS1054-03) 
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1. INTRODUCTION 



The star-formation histories of early-type galaxies in distant 
clusters have been the focus of many recent studies (e.g., v an 
Dokkum & F ranx 1996, Pender, Ziegler, & Bruzual 1996|, E l- 
hs et al. 1997, [Kelson et al. 1997[ giegler & Bender 1997|, van 
Dokkum et al. 1998 , JPahre, de Carvalho, & Djorgovski 1998 



Kelson et al. 2000c^ pegler et al. 2000[ ). While there is broad 



consensus is that E/SO galaxies in distant clusters generally 
have old stellar populations, these results depend on the ini- 
tial mass function (IMF), the cosmology, and the assumption 
that morphological evolution can be ignored. Whereas some of 
these effe cts can be taken into account using detailed mo dels 
(see, e.g., [ van Dokkum & Franx 2001 , Diaferio et al. 200C ), the 



interpretation is ultimately limited by the accuracy of the pop- 
ulation synthesis models upon which any conclusion is based. 

These assumptions can now be tested by measuring absorp- 
tion line strengths as a function of redshift (time) in early-type 
cluster galaxies. The "line strength" results can be compared to 
the evolution of colors and M/L ratios, and the correlations be- 
tween various age indicators can be compared to predictions of 
stellar population synthesis models. Together they can signif- 
icantly increase our confidence in the deduced star-formation 
histories, and in the models themselves. 

In this Letter we present first results on a new study of 
galaxy absorption line strengths to redshifts approaching unity, 
and compare the measured evolution with the evolution of the 
galaxy M/L ratios, as der i ved from the fundam ental plane (see 

and references 



present the first results using E/SO galaxies in the clusters Abell 
2256 (z=0.06), CL1358+62 (z=0.33), MS2053-04 (z=0.58), 
and MS 1054-03 (z=0.83), and discuss the resulting constraints 
on the star-formation histories, IMF, cosmological parameters, 
and the stellar population synthesis models themselves. 

2. THE DATA 

The CL1358+62, MS2053-04, and MS1054-03 spectro- 
scopic data presented in this paper were obtained at the W.M. 
Keck Observatory using the Low Resolution Imaging Spectro- 
graph ( pke et al. 1995 ). The initial sample selection for the 
deep spectroscopy was based only on magnitude and prefer- 
ence for high-resolution spectroscopic follow-up was giv en to 
spectrosc opically confirmed members ( Fisher et al. 1998 , Tran 



et al. 1999, [van Dokkum et al. 2000t [Tran et al. 2001, in prep.j ). 
Hubble Space Telescope WFPC2 images in F814W were used 
to determine morph ological classifications (Fa bricant, Franx, 
& van Dokkum 2000, van Dokkum et al. 2000, Fabricant et al. 



2001, in prep.). In this paper we limit the discussion to the E/SO 
galaxies, with more thorough discussion of the general cluster 
populations being prepared for a later paper 

The CLI 358+62 spectroscopic data and their reduction were 
published in Kelson et al. (2000b). The MS2053-04 and 
MS 1054-03 data were similarly processed. The spectroscopic 
reductions described in Kelson et al. (2000b) provided data 
of the quality required for measurement of absorption line 
strengths, with one additional step for removing atmospheric 
H2O absorption at 6800 A, 7600 A, and 8200 A. In each slit- 
mask we had included bright blue stars in order to accurately 
determine the telluric absorption in each exposure. The galaxy 
spectra presented in this paper have S/N ratios ranging from 

'Based on observations obtained at the W. M. Keck Observatory, which is operated jointly by the California Institute of Technology and the University of 
California. 

-Department of Terrestrial Magnetism, 5241 Broad Branch Rd. NW, Washington, DC 20015. Current address: OCIW, 813 Santa Barbara St, Pasadena, CA 91 101; 
kelson@ociw.edu 

'UCO/Lick Observatory, Univ. of California, Santa Cruz, CA 95064 

''Leiden Observatory, P.O. Box 9513, NL=2300 R.A. Leiden, The Netherlands 



van Dokkum et al. 1998, Kelson et al. 2000c 



therein). With this sample, we extend similar work by Ben- 
der et al. (1998) by a factor of two in redshift. Here we 
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Fig. 1.— Representative spectra for galaxies in Abell 2256 (z=0.06), CL1358+62 (z=0.33), MS2053-04 (z=0.58), and MS1054-03 (z=0.83). Tlie spectra, 
smoothed using a box width of 3 pixels, have been shifted to the restframe wavelengths, to illustrate several of the line strength bandpasses of the Lick/IDS system. 



20-100 per A (in the continuum). The resolution was typically 
o'insi ~ 40-60 km/s. The data for the local comparison sample of 
9 early-type galaxies in Abell 2256 were obtained at the KPNO 
4m and these data were processed and analyzed using the same 
procedures as were used with the high-z data in an effort to min- 
imize possible systematic errors. By doing so, we avoid any 
systematic effects between o ur distant samples an d those pub- 
lishe d in the literature (e.g., iTerlevich et al. 1999i K untschner 
2000). Figure |1| shows example spectra from each of the four 
clusters. 

Velocity dispersions were measured using the direct fitting 
method of Kelson et al. (2000b). The raw velocity disper- 
sions were corrected for aperture size to values equivalent to a 
D = 3"4 aperture at the distance of Coma using the prescription 
of J0rgensen et al. (1995). The corrections were -1-2.3%, +4.7%, 
+5.6% and H-6.0% for the galaxies in A2256, CL1358H-62, 
MS2053-04, and MS 1054-03, respectively. 

Absorption line strengths were measured using the Trager et 
al. (1998) definitions for the Lick/IDS indices. For this letter, 
we use the indices of the B aimer lines HS and H7 defined by 
Worthey & Ottavianni (1997), and will discuss the other in- 
dices in a later paper The raw indices were corrected for both 
the instrumental broadening using the Worthey & Ottavianni 
(1997) estimates for the resolution of the Lick/IDS data, and 
for Doppler broadening. 

Accurate estimation of the errors in absorption line strengths 
requi res knowledge of the variance in the data due to noise (see. 



e.g., Gonzalez 1993, Cardiel et al. 1998). In work on distant 



(i.e., highly redshifted) galaxies, the photon statistics and elec- 
tronics noise is supplemented by localized sources of noise, 
such as sharp residuals from the subtraction of the bright OH 
emission lines of the background, and residual fringing not re- 
moved by the flat-fields. Therefore, we estimate the variances 
using the residuals between the observed spectra and the high- 
resolution model SED from Vazdekis (1999), suitably broad- 
ened to match the instrumental resolution and Doppler broad- 
ening of the galaxy spectrum, which provides the lowest (as 



a function of age and [Fe/H]). We assume the square of the 
residuals between the observed and model SEDs in each index 
and continuum bandpass accurately reflects the variance due to 
noise in the data, and the errors in the absorption line strengths 
are derived by propagating these variances. Monte Carlo sim- 
ulations have verified that this method provides accurate error 
estimates. 

Because the metric apertures from which the galaxy spectra 
were extracted are systematically larger at high redshift than 
at low redshift, the B aimer line strengths needed to be cor- 
rected to a consistent aperture size. Unfortunately no published 
set of aperture corrections exist yet for the H'^a and H6a line 
strengths. However, one can estimate the corrections by in- 
tegrating the published m ean major-axis H'7 4 gradients of the 
E/SO galaxies in Fornax ( ^untschner 1999 ), weighted by an 
/•'/"^-law surface brightness profile, over a circular aperture. Uti- 
Uzing the Worthey (1994) and Vazdekis et al. (1996) model 
prediction that H5yi+H7yi = 1.9xH7/i, we find that the required 
aperture corrections scale with effective aperture diameter as 
A(H(5^+H7/i) = (1.78 ±0.16)AlogDfl,,. One can also scale the 
Mg2 aperture correction prescription of J0rgensen et al. (1995) 
by the ratio of the gradients in (HSa+Hja) and Mg2 and recover 
the same scaling. Using the above scaling of line strength with 
aperture size, we correct the line strengths to a nominal aperture 
of size D = 1"23 at the distance of CL1358H-62 (z=0.33). The 
resulting coiTections were -0.57;|]°^A, 0.17!|]{J[A, 0.20!|]{J^A, 
for A2256, MS2053-04, and MS 1054-03, respectively. 



3. EVOLUTION OF THE BALMER LINE STRENGTHS 

In Figure B we show the correlation between (HSA+H'yA) 



(Worthey & Ottaviani 1997) and velocity dispersion for E/SO 
galaxies in A2256, CL1358H-62, MS2053-04, and MS1054- 
03. For the large sample of early-type galaxies in CL1358-I-62, 
there is a well-defined correlation. This correlation spans a 
wide range of velocity dispersion and Balmer line strengths, 
from positive values (strong Balmer absorption) to negative val- 
ues (weak Balmer absorption) which occur when the mean flux 
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Fig. 2.— The (H<5a+H7a)-o- relations for early-type galaxies in A2256 (z=0.06), CL1358+62 (z=0.33), MS2053-04 (z=0.58), and MS1054-03 (z=0.83). The 
best-fitting linear relation to the CL1358+62 early-type galaxies has a remarkably low scatter of 13% in velocity dispersion. While the scatters seen in A2256 and 
MS 1054-03 are consistent with the observational errors, the observed scatter in MS2053-04 is ~50% larger than expected from measurement errors alone. The 
current dataset indicates no statistically significant differences between ellipticals and SOs. in any of the clusters. In all four panels, we show the CL1358+62 relation 
with its zero-point error (1-cr). Freezing the slope of the relation, we indicate the best-fit relations in the other two clusters using the solid lines, with the zero-point 
uncertainties indicated by dash lines. Note the mean shift to larger values of YlSA+H'yA with redshift. Increasing values of these indices are the result of stronger 
Bahner line absorption at earlier epochs. The high order Balmer indices evolve significantly with time, providing empirical calibration of the index's sensitivity to 
age. 



levels in the continuum side-bands are lower than the mean flux 
level within the index passbands. 

The form of the best-fit linear relation in CL1358H-62, shown 
with its ±l-cr zero-point uncertainties, is (HSA + H'yA) = (-9.1 ± 
1.0)log[CT/(150 km/s)] -(4.9 ±0.1). The CL1358H-62 E/SOs 
exhibit a remarkably low scatter of 0.50A, equivalent to 13% 
in velocity dispersion. When expressed in velocity disper- 
sion, this scatter is a s low as that in the fundamental plane 
( [Kelson et al. 2000c| ). We have also derived the CL1358H-62 
H7A-CT relation and have compared it to the relation for For- 
nax E/SOs (using the published data of Kuntschner 2000 ). We 
find H7yi oc (-3.6 ± 1.4)logcr at z=0.33 compared to Hja oc 
(-4.0 ± 1.5)logcr in Fornax, with similarly low scatter. The 
lack of strong evolution in the slope of this correlation sug- 
gests that stellar population ages do not vary strongly along 
the seque nce of earlv-tvpe galaxies (similar to the c onclusions 
drawn b y [Stanford, Eisenhardt, & Dickinson 1998L K elson et 
al. 2000c). 

As can be seen in Figure ^, the moderately large observa- 
tional uncertainties and/or the small sample sizes in the cur- 
rent dataset prevent an accurate measurement of the slope of 
the correlation between Balmer absorption and central velocity 
dispersion in A2256, MS2053-04 and MS 1054-03. 

Given the similarity in the slopes of the Balmer Balmer line- 
velocity dispersion relations for Fornax and CL1358H-62, we 
adopt the slope of the CL1358H-62 relation for the remainder of 
the analysis and assume that only the zero-point evolves with 



redshift. We determined the mean zero-points in the four clus- 
ters giving all points uniform weight, using the bi-weight loca- 
tion estimator (Beers, Flynn, & Gebhardt 1990). 



The zero-points for the four clusters, normalized to galax- 
ies with a dispersion of cr = 150 km/s, are -5.94 ±0.20, 
-4.86 ±0.13, -3.76 ±0.38, and -2.72 ±0.41 (in units of A). 
The A2256, MS2053-04 and MS 1054-03 early-type galaxies 
are offset in zero-point from CL1358H-62 by AiHSA + UjA) = 
-1 .08 ± 0.24 A A{HdA + H7a) = 1 . 10 ± 0.40 A and 2. 14 ± 0.43 
A, respectively. The offsets (see Figure^) are significant, 4-, 3- 
and 5-cr, and correspond to larger absorption at earlier epochs 
(i.e., younger ages), as expected. 

Figure |] shows the evolution of the (HSa + H'Ja) zero-point 
as a function of redshift. Assuming the galaxy populations in 
all three clusters can be compared directly, we show the pre- 
dictions of single burst stellar popula tion models from Worthey 
(1994; updated in Trager et al. 200C) (shaded) and Vazdekis et 
al. (1996) (dashed lines). Models with star-formation redshifts 
of 2 < z/ < 10 are shown. 

Because the model predictions can be simplified to functions 
of the form, e.g., H-ja = A-l-B&-l-(Cfa-l-D)logf -l-£'logZ-l- 
F(logZ)^, to a high degree of accuracy, we can decouple the 
time-evolution of the indices from the metallicity of a given 
model (while also rendering the adoption of Hq = 75 km/s/Mpc 
unimportant). We find lower-limits of Zf > 2.9 and Zf > 2.4 
(95% confidence limit) using the Worthey and Vazdekis mod- 
els respectively, for fi^ = 0.3, Oa = 0.7, similar to conclusions 
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Fig. 3. — The (^HA^-i-H-Y^^-rr 7prn-nnint evolution as a function of redshift, with single-burst stellar population models of Vazdekis et al. (1996) (dashed lines) 
and Worthey (1994; Trager et al. 200C ) (shaded). The model curves, shown for formation epochs of 2 < zy < 10, are insensitive to the shape of the IMF. Together 
the data and models provide lower 95% confidence limits for the mean epoch of star-formation of zf > 2.9 and z/ > 2.4 using the Worthey and Vazdekis models, 
respectively, for Q,„ = 0.3, = 0.7. 
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Dokkum et al. 1998). The predicted correlations ot single-burst stellar population models ot Vazdekis et al. (lyyo) (dashed lines) and Worthey (iyy4) (shaded) are 
also shown. The regions are defined by a range of slopes for the initial mass function 0.85 < A' < 1.85. Note: this projection of the models only assumes constant 
metallicity at fixed age and velocity dispersion, and does not depend on the assumption of co-evality. 



drawn by past work (e.g. , [van Dokkum & Franx 1996 , Bender, 
Ziegler, & Bruzual 19967Relson et al. 1997|, van Dokkum et 
al. 1998). As is always the case in this type of study, the con- 
clusions depend on the assumptions that the individual clus- 
ter samples are co-ev al, and that morphological evo lution can 
be ignored (see, e.g., [van Dokkum & Franx 2001 ). Further- 
more, we also remind the reader that early-type cluster galaxies 
may settle into their final configurations well after their last ma- 
jor epochs of star-formation leaving their morphological evo- 
lution decoupled fr om their star-formation histories (e.g., van 
Dokkum etal. 1999). 

4. COMPARISON WITH THE EVOLUTION OF M/L RATIOS 

In past work we have established the evolution of the M/Lb 
ratio for the same clusters using the fundament al plane rela- 
tion (iKelson et al. 1997i |van Dokkum et al. 1998|, Kelson et al. 
2000c) and we now compare the mean M/Lg evolution with 
that of the Balmer lines. The results are shown in Figure ^. 
Since the M/L ratios are sensitive to cosmological parameters, 
we show the results for several values of f7„, and 

The models agree very well with the observed correlation be- 
tween the evolution of the Balmer line strengths and the M/Lb 
ratios. The slope of the predicted relation depends on the slope 
of the IMF through the M/L ratio evolution (e.g., Tinsley & 
Gunn 1976). Therefore, we show the Worthey (shaded) and 
Vazdekis (dashed lines) model predictions for a range of IMF 
slopes (0.85 <x< 1.85) about the Salpeter (1955) value. For 



standard IMFs and low values of fl,„, the models are consistent 



with our data (and the results of Bender et al. 1998) 



While late-time star-bursts and/or morphological evolution 
can bias the o bserved redshift evolution of the various age in- 
dicators (e.g., van Dokkum & Franx 2001), each age-sensitive 



observable is affected similarly. Thus, comparisons like that 
in in Figure ^ provide tests of the population synthesis mod- 
els which are free from any assumption of coevality, and are 
more powerful than simple examinations of galaxy properties 
as functions of redshift alone. 

By extending our analysis to more spectral indices and galax- 
ies of later morphological type, we will better constrain the stel- 
lar and morphologi cal histories of cluster galaxies (see, e.g.. 
Kelson et al. 2000c). By doing so, we expect to quantitatively 



constrain the formation histories of distant galaxies, and better 
understand the evolution of cluster populations as a whole. 
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